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Introduction 

This report describes progress on this contract during the period 1 November 
1973 to 1 January 1974. Analysis of imagery from 11/29/72 and 1/4/73 is^entinuing. 
Digitized data from microdensitometer scans of imagery taken on 4/4/73 has been 
received and is presently being organized into a format suitable to subsequent 
analysis. Work has been initiated to locate or develop a program to read the computer 
compatible tapes for ERTS-1 imagery from 4/4/73. 

Discussion 

A majority of this reporting period was devoted to continuing the analysis of 
imagery from 11/29/72 and 1/4/74. This analysis consisted of two phases. The first 
phase involved the computation of polynomial coefficients that would yield a good 
analytical approximation to the modulus and phase constituents of the microdensitom- 
eter scan aperture optical transfer function (OTFj. The second phase analytically 
determined the required correction for irradiance falloff as a function of field 
angle (cos^^), for arbitrary scan lines, in the aerial photographs. 

In addition, the microdensitometer scans of imagery from 4/4/73 were completed 
and received. The analog output (microdensitometer chart record) has been checked 
against data log to verify order of images scanned and establish required computer 
program parameters for the data organization programs previously developed. The 
application of these programs has been initiated to organize data in a format suitable 
for analysis. 

A problem exists in the location of a program to read the computer compatible 
tapes (CCT’s) from ERTS-1 imagery from 4/4/73. This is discussed along with a 
review of the analysis phases in the following sections of this report. 
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Slit Transfer Functions 

A polynomial fit was applied to the sampled Optical Transfer Function (OTF) 
data for the microdensitoroeter scan apertures previously given in Progress Report 
#6. The resulting modulus and phase transfer functions are shown in Figures 
1-2 and supercedes all previous aperture OTF data. 

Note that the modulus is an even function of frequency so that only the even 
power coefficients of the polynomial approximations are significant . Similarly 
the phase is an odd function and thus only odd powers are used in the approximation. 
The approximating curves are therefore of the form: 
modulus M(f) = 1 + Cf^ + Ef** 

phase P(f) = Bf + Df3 

The two frequency values necessary to determine the two coefficients were chosen 
to give a good polynomial approximation over the entire frequency range of inter- 
est. In Fig. 2 two curves are shown for the 23.5 x 138 ym slit. Each curve used 
a different pair of frequency values to calculate the coefficients, giving a 
different polynomial approximation in each case. Also note the phase parts of the 
OTF's are essentially zero in the frequency range of interest. 

These curves represent the correction that must be made, in frequency space, 
to the digitized imagery data obtained from the convolution of the microdensitometer 
aperture with the photographic imagery. These aperture OTF’s will be divided into 

I 

the Fourier transform of the digitized image data to obtain the true data spectrum. 
This will then be inverse Fourier transformed to obtain the true data, as described 
in Progress Report #6. 
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Cos^ Correction 

The irradiance falloff across the image plane of an aerial camera system is 
a function of the field angle, and obeys the relationship Since, in the 

U-2 imagery scans only , a large range of field angles is covered by the scan a 
correction must be made to the exposure values. The required correction is obtained 
by multiplying these exposure values by the factor l/cos^4)(P), where P is the image 
point of interest. Since our scans are at an arbitrary angle and position within 
the format, the following process must be followed. 

Given an arbitrary scan line in the image format, as shown below, image 



Arbitrary 
scan line 


position can be described in terms of P, the distance along line PjP 2 fi^om the 

point Pi- This distance then defines r, the radius from the nadir point to the 

\ 

position P in the image format: 

1 / 

r(P) = [P^ + ri^+ 2 P (xjcose + yi sin 6)] ^ 




Then, using the geometry for aerial camera system shown below, 



the field angle, (J>, can be calculated using 


<}) = tan"^(jj 

where f is the known focal length of the camera system. 

The variation for the required correction factor, l/cos^4>, as a function of 
the image position along the microdensitoraeter scan line is shown in figures 3 
and 4 for the imagery obtained on 11/29/72 and 1/4/73. Note the significant 
difference in the value of the correction factor required between microdensitometer 
scan line 1, 2 and 3 for the imagery on 11/29/72. This is due to the fact that 
scan lines 1 and 2 were parallel and closely spaced while scan line 3 was at a 
completely different orientation and position. in the image format. This same 
argument explains the difference between the correction required for the imagery 
taken on 11/29/73 and those taken on 1/4/73, 



ERTS Computer Compatible Tapes [CCT's) 

An initial investigation into locating a program to read the ERTS CCT's has 
not been successful. Several of the ERTS investigators that use ERTS CCT’s 
were contacted during this reporting period. However, all of the programs 
utilized by these investigators were written for IBM computer systems and based 
on a 9 track tape input. The only computer system presently available for our 
use is a CDC 6400, with a 7 track tape input. 

We are presently awaiting a list of all CDC 6400 installations from Control 
Data Corporation. This will be cross referenced against a list of ERTS investi- 
gators to identify those that may be able to provide the program information 
desired. 

If our search is unsuccessful then a program will be written specifically 
for our purposes during a future reporting period. 
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Current Status of Acquired Imagery 





Aircraft 

(A/a 

ERTS- 

1 Microdensitometer 

Set 

Flight Date 

Vinten 

Scanner 

# Frames 

MSS 

scans 

1 

8/22/72 

8/23/72 

Arizona 

/ 

NA 

184 

/ 

- 

2 

11/29/72 
San Francisco 

/ 

/* 

18 

/ 

/ 

3 

1/4/73 

San Francisco 

y/** 

/* 

51 

/ 

/ 

4 

4/4/73 

San Francisco 

/ 

A 

48 


/ 

5 

6/15/73 
San Francisco 

/ 

/ 

75 

j/ *** 

6 

Requested for 

on order 

not requested 

ID not 
yet 

available 


* Scanner data not suitable for analysis because of severe geometric distortion 
arising from the lack of a gyrostablized platform on the A/C. 

** Band 001 [green) malfunction, no imagery. 

***ERTS“1 imagery of San Francisco was not acquired on 6/15/73. Therefore imagery 
from the previous [5/28/73) and following (7/3/73) cycles and from the adjacent 
6/14/73 pass have been acquired. 


Current Status of Analysis 

Calibration, matching, and scaling are in progress for the data from 11/29/72 
and 1/4/73 and all of the data is at the stage of generating OTF's. CCT's of 
ERTS-1 image # 1255-18183 [4/4/73) have been received. This data will be compared 
to the system corrected MSS photographic product to determine the amount of image 
degradation introduced in the EBR-photographic steps. 
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ERTS Frames Studied 


Data Set 

Frame #^s 

Bands 

11/29/72 

1129-18181 

4, 5, 

1/4/73 

1165-18173,-18175 

5, 6 




Figure 4: Correction 
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This report describes progress on contract number NAS 5-21849, during the period 
11/1/73 - 1/1/74. A review of the analysis of ERTS-1 and U-2 imagery is given in 
terms of correction for irradiance falloff as a function of field angle in aircraft 
images and correction of digital data by the microdensitometer optical transfer 
functions. In addition, microdensitometer scans of imagery from 4/4/73 were 
completed and received. Optical Transfer function analysis is continuing for 
imagery from 11/29/72, 1/4/73, and 4/4/73. 


